Photon-pion transition form factor: BABAR puzzle is cracked 



A. E. Dorokhov 

Joint Institute for Nuclear Research, 
Bogoliubov Laboratory of Theoretical Physics, 
141 980 Dubna, Moscow region, Russian Federation; 
Institute for Theoretical Problems of Microphysics, 
Moscow State University, RU-119899, Moscow, Russian Federation 

(Dated: April 20, 2010) 

Abstract 

Recently, the BABAR collaboration published (arXiv:0905.4778) data for the photon-pion tran- 
sition form factor i^-y-y* {Q 2 ), which are in strong contradiction to the predictions of the standard 
factorization approach to perturbative QCD. Immediately afterwards, two mechanisms were sug- 
gested (A.E. Dorokhov, TarXiv: 0905 .45771 ; A.V. Radvushkin. larXiv:0906.0323p . that logarithmically 
enhance the form factor asymptotics and therefore provide a qualitatively satisfactory description 
of the BABAR data. However, the physics of the BABAR effect was not fully clarified. In the 
present work, based on a nonperturbative approach to the QCD vacuum and on rather universal 
assumptions, we show that there exists two asymptotic regimes for the pion transition form factor. 
One regime with asymptotics F 7r7 * 7 (Q 2 ) ~ l/Q 2 corresponds to the result of the standard QCD 
factorization approach, while other violates the standard factorization and leads to asymptotic 
behavior as F^-f (Q 2 ) ~ m (Q 2 ) /Q 2 ■ Furthermore, considering specific nonlocal chiral quark 
models, we find the region of parameters, where the existing CELLO, CLEO and BABAR data for 
the pion transition form factor are successfully described. 
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I. INTRODUCTION 



In the years 1977-1981, the theory of hard exclusive processes was formulated within the 
factorization approach to perturbative quantum chromodynamics (pQCD) 1-7]. The main 
ingredients of this approach are the operator product expansion (OPE), the factorization 
theorems, and the pQCD evolution equations. In this context, the form factor for the photon- 
pion transition 7*7* — > 7r°, with both photons being spacelike (with photon virtualities 
Q\-,Q\ > 0)) was considered in jf], Q]. Since only one hadron is involved, the corresponding 
form factor F„- 7 » 7 » (Qf, Q\) has the simplest structure for the pQCD analysis among the hard 
exclusive processes. The nonperturbative information about the pion is accumulated in the 
pion distribution amplitude (DA) ip n (x) for the fraction x of the longitudinal pion momenta 
p, carried by a quark. Another simplification is, that the short-distance amplitude for the 
7*7* — > 7T° transition is, to leading order, just given by a single quark propagator. Finally, 
the photon-pion form factor is related to the axial anomaly [8j, |9j , when both photons are 
real. 

Experimentally, the easiest situation is, when one photon virtuality is small and the other 
large. Under these conditions, the form factor F„- 7 * 7 (Q 2 ,0) was measured at e + e~ colliders 



by CELLO [lOj, CLEO ll|] Collaborations (Fig. [p. In the region of large virtualities 
Q 2 » 1 GeV 2 , the pQCD factorization approach for exclusive processes predicts to leading 



order in the strong coupling constant 



)n am 



f^ D (Q 2 ,0) = ^J, (l) 

where 

J= f'dx^l (2) 
Jo x 

is the inverse moment of the pion DA, and f n = 92.4 MeV. The factor 1/Q 2 reflects the 
asymptotic property of the quark propagator connecting two quark-photon vertices (Figs. 
|2^i and [3]). The formula (JT|) is derived under the assumption, that the QCD dynamics 
at large distances (the factor Jf n ) and the QCD dynamics at small distances (the factor 
1/Q 2 ) is factorized. Moreover, under this assumption, the asymptotics is reached already 
at the typical hadronic scale of a few GeV 2 . The pion DA ip n (x), in addition, evolves 
in shape with the change of the renormalization scale (4, 6| and asymptotically equals {^J 
^As ^ = 6x (1 — x). From this follows the famous asymptotic prediction (the short-dashed 



line in Fig. [TJ) 

F^ As (Q 2 ,0) = ^. (3) 

To describe the soft nonperturbative region of Q 2 , a simple interpolation between the Q 2 — ?• 
and Q 2 — > oo limits has been proposed by Brodsky and Lepage (BL) (the dashed line in 
Fig. HD 
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Recently, the BABAR collaboration published new data (Fig. [T]) for the 77* — > 7r° 
transition form factor in the momentum transfer range from 4 to 40 GeV 2 12] . They found 
the following puzzling result: At Q 2 > 10 GeV 2 the measured form factor multiplied by the 
photon virtuality Q 2 F nl * 1 (Q 2 , 0) exceeds the predicted asymptotic limit (E} and, moreover, 
continues to grow with increasing Q 2 . This result is in strong contradiction to the predictions 
of the standard QCD factorization approach mentioned above. The BABAR data very well 
match the older data obtained by the CLEO collaboration in the smaller Q 2 region, but 
extend to a much lager Q 2 values. There is numerous literature discussing the BABAR 
effect. We refer here on, y to the first two pubficatioos fl Q appeared S ooo after the 
data were announced and it is these works that are the most relevant for the following 
consideration. In these works two scenarios were suggested, that logarithmically enhance 
the form factor asymptotics and well describe the BABAR data. 



The first scenario [l3j uses the simple constituent quark model [l5j. Within this model, 
the pion transition form factor, determined by the quark-loop (triangle) diagram with a 
momentum-independent quark mass M q , is given by 

^77*(Q 2 ,0) = -4 4^-t -{2arcsin 2 (^) + X -hi 2 h±l\ (5) 

71 W ' J 47T 2 /, ml + Q 2 2 arcsin 2 ^) 1 K 2M q J 2 (3 q -l !) K> 



where (3 q = </l + 4M 2 /Q 2 . The form factor (jSJ) has correct normalization at zero photon 



virtualities, by the axial anomaly, and has double logarithmic asymptotics ln 2 (Q 2 /M 2 )/Q 2 
at large Q 2 . This asymptotics corresponds to the case when large virtuality pass through all 
three quark propagators (Fig. [2b). In [l3] it was shown that the pion transition form factor 
calculated from (jSj) with the parameter M q = 135 MeV well reproduces the BABAR data. 
However, this model has serious shortcomings. Firstly, it has an incorrect chiral limit as 
M q — > and — > 0. Secondly, the corresponding integral for the decay constant f n within 
this quark model is divergent and thus the model should be regularized for consistency. After 
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FIG. 1. The transition form factor 7*7 — > ir°. The data are from the CELLO [10( (empty squares), 



CLEO 
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(empty triangles) and BABAR (filled circles) [12[ Collaborations. The solid line is the model 
of this work, the dashed line is the Brodsky-Lepage prediction (|4]), the short-dashed line is massless QCD 
asymptotic limit ( [3]) . 



regularization, however, the double logarithmic asymptotics is lost. Thirdly, just like in the 
Nambu-Jona-Lasinio model, it uses a local 75 vertex for the quark-pion vertex and the local 
quark propagator at all quark virtualities, in contradiction with pQCD, where there is no 
75 operator, no pion as a bound state and no constituent quark mass. It is also well known, 
that in the local quark model the distribution amplitude and distribution function of the 
pion are constants 16l-ll9|. 

Such flat (almost constant) pion DA was used in [ijj] in the context of the explanation of 
the BABAR data. The photon-pion transition form factor was calculated by using expression 
from [6] and incorporating a light-cone wave function ^ (x, k±) that has rapid falloff with 
respect to the light-front energy combination k\/x(\ — x). Within this approach, for a 
Gaussian shape of the light-cone wave function and assuming a flat pion DA <p n (x) — 1, the 
pion transition form factor is given by 



2U_ 

3Q 2 



dx 



x 



exp 



xQ< 



2(7(1 - X) 



(6) 



n 



and has logarithmically enhanced asymptotic behavior ~ In (Q /2a) /Q . In [20j it was 



demonstrated, that the descriptions of the BABAR data in the model ([6]) and in the model 
([5]) practically coincide, if o = 0.48 GeV 2 . In [14J it was also noted, that the use of such a 
wave function is numerically close to the leading-order pQCD expression for the photon-pion 
transition form factor with a modified quark propagator and a flat pion DA 

F^(Q 2 ,0) = -U^ X Q2 + M 2' ( 7 ) 

giving logarithmically enhanced asymptotics F^,(Q 2 , 0) ~ In (1 + Q 2 /M 2 ) /Q 2 . With 
M 2 = 0.6 GeV 2 the BABAR data are well fitted. Another very important feature of 
the scenario |l_4| is, that it was argumented, that there is no pQCD evolution modifications 
in the shape of flat pion DA. 

However, these approaches did not give an answer to the following serious questions. 
First of all, both expressions fl6]) and ([7]) do not describe the full form factor, but only the 
leading asymptotic part. They have incorrect normalization at Q 2 = 0, and are valid only 
at large photon virtuality Q 2 » 1 GeV 2 . Furthermore, the appearance of the parameter 
M in the asymptotic formula ((7|) is not justified. Moreover, as it was emphasized in fl^ . 
the expression ([7]) generates an infinite tower of false higher twist corrections (M 2 /xQ 2 ) n 



in contradiction with OPE. It is well known 2jJ, that there are only two terms in the OPE 



for the handbag diagram for the pion transition form factor: the twist-2 and the twist-4 
terms. Also, the relation of the parameter a in ([6]) or M in ([7]) to the fundamental QCD 
parameters and their values remained unclear. The decay constant f v is external parameter 
in this approach and is not calculated dynamically. Finally, the origin of the flat DA is not 
well understood. Most of the QCD sum rule and the instanton model calculations lead to 



the endpoint suppressed amplitudes (see, e.g. [22|, |23|). Below we show, how to generalize 
the results ©-([7]), and how to avoid the above mentioned problems with the interpretation 
of the BABAR data. 

There are several QCD based approaches to treat the nonperturbative aspects of strong in- 
teractions. They are the lattice QCD, QCD sum rules, Schwinger-Dyson approach, Nambu- 
Jona-Lasinio model, etc. In the present paper, we analyze the photon-pion transition form 
factor in the gauged nonlocal chiral quark model based on the picture of nontrivial QCD 
vacuum. The attractive feature of this model is, that it interpolates the physics at large and 
small distances. At low energy, it enjoys the spontaneous breaking of chiral symmetry, the 
generation of the dynamical quark mass, and it satisfies the basic low energy theorems. At 
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FIG. 2. The triangle diagram in momentum and a-representation notation and its different hard regimes 
for asymmetric kinematic. Hard photon is fat wavy line, real photon is thin wavy line; the fat line is the 
hard propagator, the thin line is the soft propagator: a) the standard factorization regime, b) the regime 
violating standard factorization, c) the double logarithmic regime with constant quark masses. 

energies much higher than the characteristic hadronic scale, it becomes the theory of free 
massless quarks (in chiral limit). 

The paper is organized as follows: In Sec. II, we give the basic elements of the ef- 
fective chiral quark model, the quark propagator and the quark-photon and quark-pion 
vertices. In Sec. Ill, we transform the expression for the pion transition form factor into 
the a-representation and analyze, under rather general requirements on the nonperturbative 
dynamics, the asymptotic behavior of the form factor for different kinematics. Considering 
the kinematics when one photon is virtual and other is real, we show that two possible 
behaviors of the quark-pion vertex at large quark virtualities results in two different asymp- 
totic regimes for the pion form factor. One of them corresponds to a standard factorized 
scheme with actual 1/Q 2 asymptotics. The other provides a nonstandard asymptotic regime 
leading to ~ In (Q 2 ) /Q 2 large-Q 2 behavior of the pion form factor. In Sec. IV, we specify 
two kinds of nonlocal chiral quark model implementing different asymptotic regimes and 
obtain the pion DA for various sets of parameters. In Sec. V, we are looking for the space 
of parameters that give a satisfactory fit of the CELLO, CLEO and BABAR data. Sec. VI 
contains our conclusions. 
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II. NONLOCAL CHIRAL QUARK MODEL 



Let us discuss the properties of the triangle diagram (Fig. [2]) within the effective ap- 
proach to nonperturbative QCD dynamics. To consider the asymptotics of the photon-pion 
transition form factor, we do not need to completely specify the elements of the diagram 
technique, which are, in general, model dependent, but shall restrict ourselves to rather 
general requirements. All expressions will be treated in Euclidean space appropriate for the 
nonperturbative physics. The nonperturbative quark propagator, dressed by the interaction 
with the QCD vacuum, is 

k + m (k 2 ) 



S(k) 



D(k 2 ) ' 

The main requirement to the quark propagator is, that at large quark virtualities one has 



(8) 



s(k) fe2 4°°A. 

k z 



(9) 



We assume also, that the dynamical quark mass is a function of the quark virtuality k 2 and 
normalized at zero as 

m(0) = M q , D(0) = M 2 . (10) 

At large virtualities, it drops to the current quark mass m curr faster than any power of k~ 2 
(see the discussion in [261 ]) 



m 



(k 2 )~M q exp(- (k 2 ) a ) + 



a > 0. 



(11) 



This is, first 
condensate 



because the dynamical quark mass is directly related to the nonlocal quark 
I, 0] and, secondly, the quark propagators with powerlike dynamical mass 
induce false power corrections that are in contradiction to OPE. On the other hand, the dy- 
namical quark mass (lllj) generates exponentially small corrections, invisible in the standard 
OPE. The direct instanton contributions provide a famous example of these exponential cor- 



28 



rections in the QCD sum rules approach 
is D (k 2 ) k k 2 and the typical expression is 

D (k 2 ) =k 2 + m 2 (k) . 



291 ] . The denominator in (jSJ) at large virtualities 



(12) 



It is well known (see, e.g., 30|, |31|), that the change of the quark propagator leads to a 



modification of the quark-photon vertex in order to preserve the Ward-Takahashi identity 
T M (k, q,k' = k + q) = -ie q [f - Ar M (k, q,k' = k + q)) , (13) 



where the extra term guarantees the property 

g M r M (k, q,k' = k + q) = S' 1 (k 1 ) - S' 1 (k) . (14) 

The term AT^ (q) is not uniquely defined, even within a particular model, especially its 
transverse part. The importance of the full vertex is, that the axial anomaly is reproduced 



32( | . and thus the photon-transition form factor correctly normalized. Fortunately, due 
to the fact, that AT^ is not proportional to 7 M matrix, the corresponding amplitude has 
no projection onto the leading twist operator. Thus, this term is suppressed, if a large 
photon virtuality passes through the vertex, and hence does not participate in the leading 
asymptotics of the form factor. Its leading asymptotics results exclusively from the local 
part of the photon vertex 

r^(k,q,k' = k + q) = -ie q r. (15) 
Furthermore, we need the quark-pion vertex, 

K( P ) = ^ l5 T a F(ki,e_), (i6) 

where k + and k_ are the quark and antiquark momenta. It is important to note, that the 
quark-pion vertex function F (k+Jzz) plays a similar role in our consideration as the light- 
cone wave function \1/ (x, k±_) in [lKD]. The vertex function F k 2 ,) is symmetric in the 
quark virtualities k\ and k 2 ., and rapidly decreases, when both virtualities are large. If it 
were a function of a linear combination of the quark momenta k + and then it would led 



to a growing form factor with increasing spacelike photon momenta (see for discussions |23|). 
The spontaneous breaking of chiral symmetry ensures, that the vertex function F {k 2 ., k 2 _^ is 
a functional of the dynamical mass m (k 2 ). In particular, the vertex function is normalized 
via 

F(k 2 ,k 2 ) =m(k 2 ) . (17) 

In the following, the important feature of the vertex function F (k^_, k 2 ) will be its be- 
havior in the limit, when one quark virtuality is asymptotically large and the other remains 
finite. There are two possibilities, 



Ff(k 2 + ,k 2 _) "4°°0, (18) 
and 

F^(kl,k 2 _) k -1°° g(k 2 + ). (19) 
8 



Finally, one needs the projection of the pion state onto the leading twist operator, see 
Fig. S 

T^ As (k,q,k' =* + g)=yV». (20) 

This projection is determined by the matrix element (0 \q r y flr ) 5 T a q\ n a (p)) = —i2f 7rj ps, where 
the constant f n< ps is (here w! (u) = dm (u) / du) 

r2 N c r r F(U,U) ( . , 1 . . A 

which coincides with the square of the pion decay constant f n ,ps in the so-called Pagels- 
Stokar form 33J. However note, that the physical pion decay constant, /„-, entering the 
pion vertex ffl6|) , is calculated by using the axial vertex corresponding to the conserved axial 
current (q). It turns out that the constant /^ps and the physical decay constant f n are 
not always identical. We return to this point in Sec. IV. 

Thus, we emphasize again, that in order to analyze the asymptotic behavior of the pion 
transition form factor F ni * 1 *(Ql, Q|) by inspecting the triangular diagram, one needs to 
specify only very general properties of the transition, from soft to hard regimes of the 
quark-pion-photon dynamics encoded in (jUJ), (fl5l) . (fTBl) and (I2TI) . At the same time, the 
full dynamics ([8]), f fT3|) . f fT6|) should guarantee the low energy theorems, in particular, the 
correct normalization of the form factor by the axial anomaly 

F 7r77 (0,0) = l/(47r 2 / 7r )> (22) 

and the Goldberger-Treiman relation, connecting the quark-pion coupling g qn and the dy- 
namical quark mass M q with the physical pion decay constant = Mqjg qv . 

III. ASYMPTOTICS OF PION-PHOTON TRANSITION FORM FACTOR 

The invariant amplitude for the process 7*7* — > vr° is given by 

A (7* (<?i, ei) 7* (92, e 2 ) -»■ tt° (p)) = -ie 2 e llU(XT e1e v 2 q p 1 qZFir 1 * r (-qf, -g 2 ) , (23) 

where ef are the photon polarization vectors, p 2 = m 2 , q\ = —Q\, q\ = —Q\- In the effective 
nonlocal quark-model considered above, one finds the contribution of the triangle diagram 



to the invariant amplitude 23] 



A (p 2 ; g 2 , ql) = (p 2 ; g 2 , g 2 ) + A^ (p 2 ; g 2 , g 2 ) 

9 



where the first term contains only local part of the photon vertices 

A l ° c ql ql) = ~ie 2 ^- J ^H^l "£) (24) 
■ {tr[ij 5 S(k-)e 2 S (k+ - qi)]e 1 S(k + )} + (q 1 -H- q 2 ; ei -H- e 2 )} , 

and the second term comprises the rest 

^ nonloc (P 2 ; ql ql) = -e 2 ^- J j^F{-*i, "£) (25) 
■ {tr[i7 5 S , (A;_)S' - g x ) e 1 S , (A; + )] (e 2 , Ar -q u k + - g x )) 
+tr[i7 5 S , (A;_)e25' (&+ - gi) 5(A; + )] (ei, Ar (k + - q u -q 2 , fe_))} + (gi -H- g 2 ; £i -H- e 2 ) , 

with p = q l + q 2 , q = q\ - g 2 , h± = k ± p/2. 

As we discussed above, the leading asymptotics results from the local part of the ampli- 
tude, A loc . After taking the Dirac trace and going to Euclidian metric (d A k — > id 4 k, k 2 — > —k 2 ), 
one obtains 

e 2 N c f d 4 k 2 2 m (k+) (e 12kq2 - e 12qiq2 ) - m (k 2 ) e 12qik + m (kf) e Upk 



D (k 2 + ) D (k 2 _) D (k 2 ) 

(26) 

where k 2 = (k + - gi) 2 , and e 12kq2 = e^e^fe*^, etc. 

In order to analyze the asymptotic properties of the form factor, let us transform the 



integral in ( l2~B"j) formally into the a representation (see 34j, |35|), which is one of the basic 
methods for the study of hard processes in perturbative QCD [36j , as well as in nonpertur- 

18] . Let us define for any function F of virtuality k 2 , decaying at large 



bative quark models 
virtuality as 1/k 2 or faster, its a representation (Laplace transform) 

/"OO 

F (k 2 ) = / dae~ ak2 f (a) , F (k 2 ) ~ / (a) , (27) 
Jo 

where F (k 2 ) is the image of the original / (a). The important asymptotic property of this 
representation is, that the large power-like k 2 behavior of F (k 2 ) is given by derivatives of 
the original g (a) at a = 

Thus, the large k 2 asymptotics of the image F (k 2 ) is related to the small a behavior of the 
original / (a) . 
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Let us introduce the following notations 
1 , , , m(k 2 ^ 



D(k 2 ) 
F(k 2 + ,k 2 _ 



d (a) 



D(k 2 
G(a,(3), 



~ dm (a) 



rn 



(hi) F(k 2 +J k 2 



G m fl (a, 0) 



(29) 
(30) 



D(k 2 + )D(k 2 _) D(k\)D(k 2 _) 

where in the second line the double a representation is implied. Because of the properties 
© and dnj one has 

d(0) = l,d' (0) = 0,d"(0) = 0,..., 
d m (0) = 0,tT m (0) = 0,... 



With this notation, using the standard technique of the a representation ( 34|, [35|), the 
momentum integral in (|26|) is transformed into the following expression for the form factor 



f*v {p 2 ;QlQl) 



p -i[-^ 2 +7(qQf+Wl)] 



[d (7) (aG m>0 (a, /?) + /3G , m (a, /?)) + 7dm (7) C (a, 



(31) 



where A = a + (3 + 7 and J d (a/3j) ... = / °° da J °° d(3 J °° d7. 



A. Symmetric kinematics 

Let us first consider the symmetric kinematics Q\ = Q\ = Q 2 . Then one has 

i^ 7 V <?> Q 2 )=^fj ^^e-U-^(^} (32) 
• [d (7) (aG m ,o (a, /3) + /3G , m (a, /?)) + 7d m (7) G (a, /?)] . 

Large Q 2 behavior of F^?* * (p 2 ; Q 2 , Q 2 ) corresponds to either small 7, small a + 0, or to 
large A. It is easy to check, that the leading asymptotics is ensured by small 7 and thus 
A — > a + (3. The term with factor 7d m (7) provides only exponentially small corrections 
and does not contribute to the leading asymptotics. In this way, in (132]) . the integral over 7 
(small distances) and the integral over a, /3 (large distances) is factorized. The integral over 
7, using (1271) . transforms the original d (7) back to momentum space 1/D (Q 2 ) 

(p 2 ; Q 2 > Q 2 ) Q =°° w UD l m J ^f e ^ p2 (a ' ® + ^ (a ' M ■ 

(33) 
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For the quark propagator, one obtains in this limit 1/D(Q 2 ) — > 1/Q 2 plus exponentially 
small corrections, due to the properties ([9]) and ffTTj) . It turns out, that the integral in f )33|) 
is the a representation of the pion decay constant (12~T1) 



if 



d(a(3) ji&2 



lps >* 4n 2 J (« + /3 ) 3 
Thus one obtains the asymptotic formula 

i><L, (o ; q 2 ,q 2 ) 



2 1 /PS,7T 



(34) 



(35) 



3Q 2 U ' 

for the form factor in symmetric kinematics, which for models, where fps,n = fn, reproduces 
the Brodsky-Lepage factorization result [7|. 

(WWW 7 * (9l) 




■WWW 
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FIG. 3. At large photon virtualities in the leading twist, in the srandara facxorization regime, the amplitude 
is factorized into the soft pion matrix element and the hard coefficient function. 



In order to define the pion DA, we carry out a change of variables in 



a — > xL, (3 — > xL, 



with x — (1 — x) , then 



W (X) 



N n f°° dL ^ X - Lp2 



with 



e xxL P (xQ xL ^ _j_ a;G! (xL, xL)) 



dxip n (x) — 1. 



(36) 



(37) 



In the momentum representation and using the chiral limit p 2 = 0, the result (1371) for the 
leading twist DA is [23] 



iV c f°° dX f°° , F(u + zAx,m-2Ax) . . . , 

¥M^J = 73^72 — / 77Z / du-p^j , - _ N n ^ — . x _ x [u + ^Ax) + xm [u — ?Ax)j . 



47r2 /p S n J -00 2tt 7 """D (u - i\x) D{u + iXx) 



(38) 
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This result is also in agreement with earlier calculations made in the instanton model 



under some simplified assumptions 19|, 37 ^39(] . The arguments in the integrand have the 
simple meaning of the transverse u = k\ and longitudinal parts of the quark (antiquark) 
virtualities. 

For the pion vertex with the property (FI5|) . the pion DA vanishes at the endpoints 

V l (x = 0) = ipf (x = l) = 0, 



while for the second type of the pion vertex (fT9|) . one has instead 

Nr. 



V U J (x = 0) = V U J (x = 1) 



™ du m{u)g{u) 



D(u) 



(39) 



The pion DA ip n (x) in (|37|) is the leading twist- 2 DA, defined as a gauge- invariant matrix 
element of the nonlocal operator 







d (z) 7^7 5 P exp / dz^A^ (z) u {-z 



vr+ (p) ) = iV2f™ Pfl / dxe^^if^x) 



(40) 

with the Dirac structure 7^75 between the pion and vacuum states, z a light-like four- vector 
[z 1 = 0), and the gluon field (z). 

Thus, in symmetric kinematics, the standard factorization is not violated and the OPE 
is modified only by exponentially small terms. 



B. Asymmetric kinematics I 

Let us now consider the asymmetric kinematics Q\ = Q 2 , Q\ = 0. Then one has 

N c f djaM) 1 [_ Q;3p 2 +7aQ 2] 



(A Q\ 0) = ^j- j ^R e -±[-^^\ (41) 

• [d (7) (aG mfi (a, 0) + /3G ,m (a, 0)) + ld m (7) G (a, (3)} . 

For simplicity in the following we shell consider the chiral limit, m curr = 0,p 2 = 0. 

Let us first consider the model with the quark-pion vertex possessing the property ( fl8j) . 
In this case, the regime of small a does not lead to the leading asymptotic terms because of 
property G (a, 0) — > as a — > 0. The leading large Q 2 behavior corresponds to small 7, i.e. 
A — > a + (3, (Fig. [2^), as for symmetric kinematics, 

^ C 4(0;Q 2 ,0) Q ^°° (42) 
vn z U J ( a + /3) 
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This asymptotic term corresponds to the standard factorization contribution (Fig. [3]) and 

the integral over 7 again can be transformed back to the momentum space 

z>c,i / n r,2 n \ Q 2 ^°° N c f d(a/3) aG mj0 (a,f3) + /3G , m (a,f3) 
b wtl , OJ - 



After change of variables (1361) . we arrive at the representation 

^(O^O^f^jfd^rfM, (43) 

where <p n (x) is defined in fl37|) . Because in the considered case (p n (x) vanishes at the 
endpoints the actual asymptotics is 

F^(0-Q\0) Q2 ^±l f -^P (44) 

in agreement with (Q]), where J? = ^-(p£ (x) is given in the momentum space representa- 
tion as [23] 

If N ' r A u (\ Ft (u,yu)m(yu) 



As we have already noted in Introduction the asymptotic behavior (1441) is not seen in 
the BABAR data. Nevertheless, even for the case considered, in principle, it is possible to 
simulate in some wide preasymptotic kinematical region a logarithmically enhanced behavior 
of the form factor. This happens if one assumes that the pion DA entering (]4"3l) is almost 
flat (x) ~ 1, i.e. it is close to a constant everywhere except small vicinity near endpoints. 
Then, in order to regularize the integral for J* in the infrared region, one needs to keep the 
exponentially small terms in ( 143|) . 

To this end, let us analyze the asymptotic behavior of the integral 

jL = Q2 [ dX DW)> (46) 
corresponding to a flat pion DA, for some popular models of the nonperturbative quark 
propagator. Firstly, we consider the quark propagator 

1 1 -exp(-/c 2 /A 2 ) 



D (k 2 ) k 2 

with the property of analytical confinement 0, I^H- I n quark models, where this propagator 

is used, the parameter A has the meaning of a dynamical quark mass j^, A = M qj with 

typical values of M q = 200 - 300 MeV. Inserting (jSD into (USD one obtains 

f 1 l-exp(-xQ 2 /M 2 ) , s 

J L ac = I dx V q) (48) 
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with the leading asymptotic behavior 

J ac Q2 =°° In (Q 2 /M 2 ) + 7b, (19) 



where 7^ is the Euler-Mascheroni constant. Both expressions (I48p and (|49|) are very close to 
the result ([2]) obtained in [lj| . The difference is, that in the expression ( I48p the extra factor 
(1 — x)" 1 in the exponent is absent, and more important the parameter in the exponent in 
( 148]) has clear physical sense as a dynamical quark mass squared. 

Secondly, let us take the propagator of the general form given in ([8]) 

J ^ Q2 [ ix ^ + lwY (50) 

Then one obtains the asymptotic behavior 

Jo (u + m 2 (u)) [u + M 2 ) 

Again, this is similar to (JZJ) obtained in [l^], but with important differences. In fact, (JTj) is a 
purely asymptotic formula and it is not allowed to keep the parameter M in the asymptotic 
quark propagator. The expression ([50]) is valid for all Q 2 and provides the leading asymp- 
totics for the flat DA (l5Tj) . It has correct large Q 2 behavior for the quark propagator, 1/Q 2 , 
and does not contain false power corrections. 

C. Asymmetric kinematics II 

Now, let us consider the model with the quark-pion vertex possessing the property ( TT9l) . 
It is convenient to rearrange the terms in the pion form factor in the following way 

+ aG mfi {a, 0) d(j) + p [G , m {a, p) - r m 
+ 7 G (a, 0) d m (7) + (3r m [d (7) - 1] 
+f3 [tZ( 7 )-l] \PG 0>m (a,P)-r m (/?)]}, 

where we introduce notations for the originals 

g(k 2 ) m{k 2 )g{k 2 ) 
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The term in the fourth line of (1521) vanishes as 7 — > or a — > 0, and thus does not 
participate in the leading asymptotics. The terms in the second line vanish as a — > 0, but 
remains finite as 7 — > 0. In this case, the 1/Q 2 asymptotics is due to hard quark propagator 
connecting two photon vertices and the coefficient reflects the soft properties of the pion (Fig. 
[2k and Fig. 3). For the terms in the third line one has an opposite situation, they vanish 
with 7, but finite as a — > 0. Thus, the 1/Q 2 asymptotics is due to hard quark propagator 
connecting pion and hard photon vertices, while the coefficient correlates soft properties of 
the pion and photon (Fig. [2b and Fig. 4). The term in the first line of (1521) provides the 
asymptotics ~ ln(Q 2 ) /Q 2 . This asymptotics corresponds to a combined soft-hard regime 
when one parameter (i.e. a) vanishes, while the other (7) goes to infinity 

After standard manipulations with the integrals one obtains the following large-Q 2 asymp- 



totic behavior transformed to the momentum representation 



F^(0 ]Q 2 ,0) Q2 -^ 1 V 



Q 2 6^ 

00 1 rl 



00 m{u)g{u). (Q 2 



du _ y J V y ln ^- ) +A 
D (u) \ u 



(53) 



A = i°° du Dj^i d yp^{ uFUf ^y u )-[ u + 2m2 ( u )]9(y^)}- (54) 

The coefficient of the logarithmic term in ([531) is clearly related to the fact that the pion 
DA for the case considered does not vanish at the endpoints and proportional to the value 
of the pion DA at these points, see ([39]) . When the function g (u) = 0, we reproduce 
the asymptotics (f44j) and (145!) corresponding to the quark pion vertex with property (|T8l) . 
The variable u in the integral ( 1531) may be considered as the square of quark transverse 
momentum in the pion, u ~ k\. The asymptotic expression (1531) generalizes the asymptotic 
formula ([1]) for the case when the standard factorization is violated. 

IV. THE INSTANTON AND CHIRAL MODELS 

In the previous section we considered the asymptotic behavior of the pion transition 
form factor given in f[2"3"|) - ([2"5'|) . In order to calculate this form factor in the whole kinematic 
region and compare with available experimental data, we should further specify our model 
assumptions. Let us introduce the momentum-dependent dynamical quark mass entering 
the propagator ([8]) as (we consider the chiral limit m curr = 0) 

mjk 2 ) = M g f 2 (k 2 ) (55) 
1 See for classification of different regimes 36| . 
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7*(gi) 



7(92 



7* (9i) 

IWWW 




FIG. 4. The factorization of the amplitude into the soft part relate' 
wave function and the hard part of the quark propagator. 

and take the profile function / (k 2 ) in a Gaussian form 

/ (k 2 ) = exp (-AA; 2 ) 



projection of the photon 



(56) 



Thus, the model contains two parameters, the dynamical quark mass M q and the non-locality 
parameter A. 

Next, we need to specify the nonlocal part of the vector vertex that does not participate 
in the leading asymptotics, but is very important in implementing the low energy theorems. 
The nonlocal part of the vector vertex in ( Tl"3l) is taken of the form 30] 



at, /, ,/ , x /, r/s m(k' 2 ) — m(k 2 ) 
AT, (k, q ,k' = k + q ) = (k + k ) fM \! 2 _ k2 K 



(57) 



Further, we will consider two kinds of quark-pion vertex (jl6j) . the first given by 



and the second by 



Fi (k\, k 2 _) = M q f (k 2 + ) f (k 2 _) 



F x (k 2 + ,e) = ^M q [f 2 (k 2 + ) + f 2 (k 2 .)} 



(58) 



(59) 



The first one is motivated by the instanton picture of QCD vacuum [24| and the second by 
the nonlocal chiral quark model advertised in [25]. We shall in the further discussion refer 
to vertex function ( 1581) . which has the k 2 — > 00 behavior (TIBl) . as the instanton model, and 
to the other choice (159]) . corresponding to k 2 — > 00 behavior (I19p . as the chiral model. 

The important requirement, that correlates the parameters of the models, is to fit the pion 
decay constant For the instanton based model this constant is given by the expression 
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found in 241 ] 



/dp,tt 



4tt 2 



m[u) 



d u u r^o ( m (m) — to' (m) + it (m' (w)) ) , 
(it) 



(60) 



and for the chiral model 



25| the expression for f n coincides with the Pagels-Stokar form 



( 12 ip . Within the nonlocal chiral model approach there is a difference between the vertex 
corresponding to the conserved axial current, 



Tl (k, q,k' = k + q) = [7 V - Arj (k, q,k' = k + q)] 



(61) 



and the local vertex (I20p . corresponding to the leading twist operator. The total axial vertex 
(q) ensures the axial Ward-Takahashi identity and the Goldberger-Treiman relation. The 



non 



31 



oca. 



45 



part of the axial vertex, that leads to (j2ip is given in 44j and to fl60l) is given in 



46]. 



Fig. shows the parameter space where the pion decay constant is fixed by its value 
taken in the chiral limit f n = 85 MeV ^7]. 



2,5 



2,0 



„ 1,5 

> 
O 

O 

<: 1,0 




0,40 



M q (GeV) 

FIG. 5. The correlation between dynamical quark mass M q and the nonlocal parameter A that fit the pion 
decay constant in chiral limit f n = 85 MeV. The solid line is for the chiral model fps,-w = fm an d t ne 
dashed line is for the instanton model Jdp.it = fir- 
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For the instanton model (1581) . the pion DA (|37|) is reduced to 
]y r°° d L 

& ( x ) = , 2 / 2 — M 9 / ~T~ ( xo ™ °" + ( xL ) a ™ (xL)) , (62) 
rf(x = 0) = 0, / dx<&(x) = l. 



For the chiral model (1551) . one obtains the pion DA 

AT Z -00 m 2 (til Z -1 
(* = 0) = — 2 c - 2 / rfn— ^ / rf^x (x) = j (64) 



4vr 2 / 2 7 D (u) 
which is not vanishing at the endpoints x = and x — 1. 

In above expressions, we used the following notations for the correspondence between 
momentum and a-representation (in addition to definitions ( 1291) . (150]) ) 

m 2 (A; 2 " 1 



D(ife 



2^ 



d m2 (a) 



/(fc 2 ) , . m(fc 2 )/(fc 2 ) 



D (/c 2 ) w ' D(A; 2 
The explicit form of the functions in a representation in the case of the model defined by 
( 1551) and ( 1561) is given in Appendix. 

In Fig. |6] the different shapes of the pion DA are shown as they are calculated within the 
instanton and chiral models for the values of the dynamical quark mass M q = 300 MeV and 
M q = 125 MeV. The parameter A is defined to fit the pion decay constant in chiral limit 
f n = 85 MeV. For smaller M q the pion DA is close to a flat shape. For larger M q it is more 
sensitive to the nonlocal part of the photon vertex and, in case of the instanton model, it is 
strongly suppressed in the vicinity of endpoints. 

In Fig. [7] the prediction for the pion transition form factor in symmetric kinematics 
calculated from (j2"3"|) - (l2"5|) is presented. The explicit expression for the instanton model is 

(0; Q\ Q 2 ) = ^J ^e-i^V (fi) (65) 

■ [2aa m (a) d (7) + 7a (a) d m (7)] , 

and for the chiral model is 

F%*. (0; Q 2 , Q 2 ) = ^J-J d ^ 1 *d (7) [e-^^d m2 (a) d (p) (66) 



+ (Vi7(* + /3)Q 2 + e -^0? + 7)Q 2 ) dm (q) dm {(3) 
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d) 




FIG. 6. Pion distribution amplitude for the instanton model with parameters a) M q = 125 MeV, A = 0.016 
GeV" 2 and b) M q = 300 MeV, A = 1.3 GeV~ 2 ; and chiral model with parameters c) M q = 125 MeV, 
A = 0.0098 GeV" 2 and d) M q = 300 MeV, A = 0.639 GeV~ 2 . 

As it is seen from FigJTl the qualitative behavior of the pion transition form factor for fixed 
quark mass is similar for the two different models. For M q = 300 MeV, the combination 
Q 2 F 7T1 * 1 * rapidly turns into the asymptotic regime as expected in the standard factorization 
scheme. The asymptotic limits are different for the two models, 2f n /3 for the chiral model 
and 2/p S7r /3/ T for the instanton model. However, for smaller masses the effect of vertex 
non- localities is diminished, in particular fop,n ~ /ps> for the instanton model. One sees 
from Fig. [3, that for M q = 125 MeV the behavior of the form factors is similar for both 
models. 

V. THE BABAR DATA WITHIN THE INSTANTON AND CHIRAL MODELS 

Let us consider the model predictions for the pion transition form factor in the asymmetric 
kinematics (qf = Q 2 ,q% = 0) calculated from ( 123]) - (T25|) in the region, where experimental 
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data exist. The explicit expression for the instanton model is 



(0;Q 2 ,0) 



N c M q r d (gpj) _^ Q 2 
~6^f n J A3 6 

{(aa m (a) a (f3) + (3a (a) cr m (/?)) d (7) + 7a (a) ex (/3) d m (7)] , 



(67) 



and for the chiral model is 



F#*(0;Q 2 ,0) 




j dj^ e -^ {idm{a)mdm{i) 



(68) 



+d (7) [arf m2 (a) d(f3) + d (a) /?d m2 (/?) + (2a + 0) d m (a) d m (/?)]} . 



In Fig. |H1 we show the predictions for different values of M g .For a quark mass M q = 300 
MeV the model dependence is very strong and the theoretical curves are very far from the 
experimental points. The chiral model overshoots the data, while the instanton model, 
in correspondence with the standard factorization scenario, shows the asymptotic 1/Q 2 
behavior very early, already at Q 2 ~ 1 GeV 2 . It is clearly seen, that in order to describe 
the BABAR data, one has to take the dynamical quark mass M q 125 MeV. Then both 
models have an qualitatively good description, with some preference to the chiral model. In 
Figs. [9^ and ITUa we show that the parameter space that describes the data up to 40 GeV 2 is 
rather narrow. For the chiral model it is M q pa 125 ± 10 MeV, and for the instanton model it 
is M q pa 130 ± 5 MeV. Thus in this region the instanton model simulate the logarithmically 
enhanced behavior due to rather flat pion DA. However, the further behavior of the form 
factor is rather different for different models as it is seen in Figs. |9b and [TOb . where the 
kinematical region up to 100 GeV 2 is shown. The instanton model finally reach its actual 
asymptotic 1/Q 2 that follows from f|44|) and ( I45p with the asymptotic coefficient given by 



For the chiral model the logarithmic growth continues for all Q 2 with the asymptotics fol- 
lowing from ( |53l) 




(69) 



F^(0;Q 2 ,0) 



Q 2 ->oo 



1 N c 
Q 2 12vr 2 /, 



00 



du 



m 2 (u) 




(70) 



7T 
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Let us make few comments. First of all, the form factor and its asymptotics are rather 
different at lower Q 2 . From Fig. [11] it is seen that the asymptotic curve conjugates the 
calculated curve in the region of order of 100 GeV 2 . Secondly, the fact, that the quark mass 
leading to a satisfactory fit of the data is quite small, is not fully unexpected. There are 
not many quantities that are very sensitive to the dynamical quark mass. The precisely 
known contribution of the hadronic vacuum to the anomalous magnetic moment of muon, 
q — 2, is infrared sensitive and demands low values for the quark mass, M q « 200 MeV 
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50] . Finally, remember also, that understanding the asymptotics of the pion transition 



form factor is important for selection of realistic nonperturbative models, used to estimate 
the hadronic contribution of the light-by-light process to g — 2 
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52|. 



VI. CONCLUSIONS 

The present paper is devoted to the so-called BABAR puzzle. New very precise data 
were obtained by the BABAR collaboration for the photon-pion transition form factor in 
very wide kinematical region up to large photon virtualities Q 2 ~ 40 GeV 2 [12j|. The data 
overshoot the asymptotic limit for Q 2 F 7T11 * (Q 2 ) predicted by Brodsky and Lepage [7|, and 
have a tendency to grow further. Both facts are in strong contradiction with the standard 
QCD factorization approach, which constitutes the BABAR puzzle. 

The main problem is the unstopped growth of the new data points for Q 2 F. K11 * (Q 2 ) that is 
inconsistent with the predicted Q 2 F nll * (Q 2 ) — > constant, following from simple asymptotic 
properties of the massless quark propagator. The key point, to solve this problem, is to 
consider the properties of the pion vertex function F{k\ 1 k 2 ) which is the analog of the light- 
cone pion wave function. There are two possibilities for the momentum dependence of the 
pion vertex function. In the limit, when one quark virtuality, k 2 , goes to infinity, and the 
other, k 2 , remains finite, the vertex function may not necessarily tend to zero. When it goes 
to zero, the pion DA (p^(x), which is a functional of the pion vertex function, is zero at the 
endpoints, ip n (0) = ^(1) = 0, with either strong or weak suppression in the neighborhood of 
the endpoints x = and x — 1. For the situation of strong suppression, the asymptotic l/Q 2 
behavior of the pion form factor in asymmetric kinematics (Qf = Q 2 ,Q\ = 0) is developed 
very early, in contradiction with the BABAR data. For weak suppression (resembling a flat 
distribution amplitude of the pion), the asymptotic l/Q 2 behavior is developed quite late, 
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and can give a reasonable description of the data in the BABAR region with a \nQ 2 /Q 2 
behavior in this region. For the other case of non-vanishing pion vertex function in the 
above limit, the pion DA ip n (x) is not zero at the endpoints, and therefore the asymptotic 
In Q 2 /Q 2 behavior persists over the whole range, in particular in the BABAR region. 

In order to fit the available data on the photon-pion transition form factor from CELLO, 
CLEO and BABAR, we have analyzed the parameter space of two examples of nonperturba- 
tive models, motivated by the instanton 24| and the chiral 25[ models, characterized by the 
two parameters, dynamical quark mass M q and the parameter of non-locality A. The main 
conclusion is, that the fit to the data requires a quite small dynamical quark mass M q w 125 
MeV with rather small uncertainty. As a consequence, the parameter of non-locality, that 
fits the pion decay constant f n , is very small, A ~ 0.01 GeV~ 2 . Thus, one has an almost local 
quark model with very flat regulators in momentum space, that considerably diminishes the 
difference between the nonperturbative models considered in this work. In this respect, this 
situation resembles the fit by a simple local quark model, made in [3] with M q = 135 MeV. 
On the other hand, in [l4] only the leading asymptotics were used and large mass parameter 
of order of 1 GeV were required, to fit the BABAR data. 

Finally we would like to point out, that in the present work, we did not consider QCD 
evolution. In [ijj], it was argumented, that the flat pion DA corresponds to a very small 
momentum scale, and hence QCD evolution is frozen. Our calculations support this point 
of view. In particular, our choice of parameters fitting the BABAR data leads to quite large 
values of the quark condensate, also corresponding to a very low normalization point. 

Concluding we may say, that the BABAR data being unique in their accuracy and cov- 
ering a very wide kinematical range, are consistent with considerations based on nonpertur- 
bative QCD dynamics and may indicate specific properties of the pion wave function. 
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VIII. APPENDIX 

Here, the explicit expressions of the functions in a representation for the Gaussian model 
defined by (153]) and (I5B]1 are given 

oo (_-t\n 

d (a) = i + ^-f- \ M l ( a ~ 4Ari )] n (" - 4An ) > 

n=l ^" 
co (_-,\n 

d m (a) = M q J2 ^"7^ IK (a - A (2 + 4n))] n (a - A (2 + An)) , 

n=0 U ' 
co /_-,\n 

rf m2 («) = M g 2 J]^-^[M g 2 (a-A(4 + 4n))] n 0(a-A(4 + 4n)), 

n=0 

a (a) = J2^~ [M ( 2 (a-A(l + 4n))] n 6(«-A(l + 4n)) ) 

n=0 n ' 

co (_-,\n 

a m (a) = M q 2J [Mj (a - A (3 + 4n))] n 6 (a - A (3 + 4n)) . 

n=0 n ' 
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Q 2 GeV 2 



FIG. 7. Photon-pion transition form factor in symmetric kinematics for the instanton model with parameters 
M q = 125 MeV, A = 0.016 GcV~ 2 (short pointed line), M q = 300 MeV, A = 1.3 GeV" 2 (dash-dotted line); 
and chiral model with parameters M q = 125 MeV, A = 0.0098 GeV~ 2 (solid line) and M q = 300 MeV, 
A = 0.639 GeV -2 (dashed line). The straight dotted line is asymptotic limit 2f w /3. 
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Q 2 (GeV 2 ) 

FIG. 8. Photon-pion transition form factor in asymmetric kinematics for the instanton model with param- 
eters M q = 125 MeV, A = 0.016 GcV~ 2 (short pointed line), M q = 300 MeV, A = 1.3 GeV~ 2 (dash-dotted 
line); and chiral model with parameters M q = 125 MeV, A = 0.0098 GeV~ 2 (solid line) and M q = 300 MeV, 
A = 0.639 GeV -2 (dashed line). The straight dotted line is asymptotic limit 2f w . 
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FIG. 9. Photon-pion transition form factor in asymmetric kinematics for the chiral model with parameters 
M q = 125 MeV, A = 0.0098 GeV~ 2 (solid line), M q = 135 McV, A = 0.0203 GcV~ 2 (dashed line), M q = 115 
MeV, A = 0.0038 GeV~ 2 (dotted line). The straight dotted line is asymptotic limit 2f n . 
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FIG. 10. Photon-pion transition form factor in asymmetric kinematics for the instanton model with pa- 
rameters M q = 115 MeV, A = 0.0038 GeV~ 2 (dotted line), M q = 125 MeV, A = 0.0098 GeV~ 2 (solid line), 
M q = 135 MeV, A = 0.0203 GeV~ 2 (dashed line), M q = 150 MeV, A = 0.077 GeV~ 2 (dash-dot-dottcd line). 
The straight dotted line is asymptotic limit 2f n . 
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